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(54) Linear space-time multiuser detector 

(57) Linear space-time multiuser detectors are 
described that combine array processing, rake detec- 
tion, and multiuser detection and that satisfy a minimum 
mean-squared en'or criteria. Both embodiments can be 
implemented adaptively so as to account for unknown 
sources of interference. The adaptive implementation of 
the first embodiment of a linear space-time detector 
requires explicit estimates of the channel coefficients. 
On the other hand, the adaptive implementation of ttie 
second embodiment does not require these estimates 
since they can be adaptively obtained. Both of these 
detectors provide significant performance advantages 
over single-user space-time detectors and over multi- 
user detectors which do not account for CDMA interfer- 
ence from unknown sources. 
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Description 

Field of the Invention 

5 [0001] This invention relates generally to direct sequence (DS) code division muitiplex access (CDMA) systems sut>- 
jected to multipath bding and muKq^le access interference (MAI), and more particularly to signal detection techniques 
for use in such systems. 

Badcground of the Invention 

10 

[0002] In wireless CDMA systems such as those proposed for next generation mobile wireless standards, the major 
impediments at the physical layer are multipath bding and multiple access interference caused by co-channel users 
which are not orthogonal to the desired user. The rake receiver, described in R. Price and P. E. Green. "A Communica- 
tion Technique for Multipath Channels," Proceedings of the IRE, Vol. 46, pp. 555-570, March 1958, attempts to con^at 
15 multipath fading by coherently combining resolvable nujltipath replicas of the desired signal. Multiuser detection 
desaibed in S. Verdu, "Multiuser Detection,'' Cambridge University Press. New York, 1998, addresses the problem of 
MAI by actively accounting for its presence when detecting the desired user. 

[0003] More recently, there has been growing interest in using array processing for further improving receiver perform- 
ance. These techniques have focused on using multiple antennas at the base station receiver to provide antenna gain 

20 and/or diversity gain and allow the possibility of spatial processing. By conrtfjining these space-domain techniques with 
time domain techniques like rake detection and multiuser detection, the resulting space-time detectors show promise 
of improving the capacity of CDMA systems as compared to traditional time-domain-only detectors. See, for exanple. 
A. Paiilraj and C. Papadias, "Space-Time Processing for Wireless Communications,' IEEE Signal Processing Maga- 
zine, Vol. 14, No. 6, pp. 49-83, Nov. 1997. The first generation of space-time CDMA detectors used anay processing 

25 with either rake detection or multi-user detection. See, respectively, A. Naguib and A. Paulraj. "Performance of Wireless 
CDMA with M-ary Orthogonal Modulation and Cell Site Antenna Anrays," IEEE Journal on Selected Areas in Commu- 
nications, Vol. 14, No 9, pp. 1770-1783, Dec. 1996 or S. Miller and S. Schwartz, "Integrated Spatial-Temporal Detectors 
for Asynchronous Gaussian Multiple-Access Channels," IEEE Transactions on Communications, Vol. 43, No. 2/3/4, pp. 
396-411. Feb.yMar7Apr. 1995. Later space-time CDMA detectors confined all three processing techniques. See: H. 

30 Huang; S. Schwartz, S. Verdi!i. "Combined Multipath and Spatial Resolution for Multiuser Detection: Potentials and 
Problenis." Proceedings of the IEEE International Symposium on Information Theory, p. 380, 1995; or M. Nagatauka 
and R Kohno, "A Spatially and Temporally Optimal Multi-User Receiver Using an Array Antenna for DS/CDMA" lEICE 
Transactions on Communications, Vol. E78-B, No. 11. pp. 1489-1497, Nov. 1995. 

[0004] While the foregoing systems operate satisfactorily, improvements can be made, in particular to tiie space-time 
35 detectors which combine all tiiree processing techniques. While tiie detector in ttie Nagatsuka and Kbhno paper is opti- 
mum in the maximum likelihood sense, its computational complexity is exponential with respect to tiie number of users. 
Hence it is too complex to implement for practical systems. In the paper by Huang. Schwartz and Verdt), a tradeoff 
between performance and complexity is made, but this detector was not implemented adaptively since it used a zero- 
forcing criteria. Adaptive implmentations allow receivers to account for unknown sources of interference ttius improv- 
40 ing the detector performance and Increasing the system capacity. For example, a base station receiver could account 
for interference from adjacent cells or from an embedded microcell, while a handset receiver oouM account for interfer- 
ence from signals it is not explidtiy demodulating. 

Summary of the Invention 

45 

[0005] This invention detects DS-CDMA signals utilizing a rake receiver, array processing, and multiuser detection. 
When combined with array processing, the rake receiver is often called a space-time rake receiver. The invention uses 
a minimum mean-squared enror (MMSE) criterion in ttie multiuser detector. This criterion allows for relatively simple 
implementation In the form of a linear detector and also allows for adaptive implementation. Adaptive implementations 

50 are useful in practical situations where there is limited knowledge of ttie various received signals. As discussed above, 
botti uplink and downlink capacity can be improved using adaptive detectors which account for unknown interference. 
[0006] For pulse anrplitude nxxf ulated (PAM) data signals, two embodiments of ttie invention provide options for trad- 
ing off between performance and adaptive implementation complexity. The first eni^odiment perfbmns better if tiiere is 
perfect knowledge of ttie users' signal parameters. On the ottier hand, it requires more explicit channel information for 

55 adaptive implementation. Suppose the receiver consists of 1 antennas and the received signal consists of K DS- 
CDMA data signals, each witti L delayedMeighted multipath replicas. At each antenna, a bank of filters is matched to 
the KL spreading codes with their multipath timing delays. The filter outputs are weighted according to the complex con- 
jugate of ttie estimated channel (muttipatti and array) parameters and combined to form a K-vector where each com- 



EP0 978951 A2 



ponent corresponds to one of the K codes. Hie real part of each component is taken. Estimates of the multipath delays 
and channel parameters may be okstained. for example, from a training or pilot signal. The conventional space-time rake 
receiver would pass each vector component to a decision device for estimating the corresponding user's PAM data 
symbol. However, unless the users* signals are orthogonal in the space-time domain, the vector components are con- 

5 taminated with multiaccess interference from the other users. This invention proposes the use of a linear conibiner prior 
to the decision devices to suppress the multiaccess interference. This linear combiner is represented by a real K-by-K 
matrix Wj^ which multiplies the real fC-vector. The matrix minimizes the mean squared error between its product with 
the real K-vector and the /C-vector of data symbols. Because the matrix uses the minimum mean-squared 
error(MMSE)criterion. well-known adaptive aigorittims can be used to adaptively obtain it Each component of the final 

10 K-vector output con^esponds to one of ttie /C codes and is passed to either a decision slicer or a decoder for further 
processing. 

[0007] The second en^odiment has marginally inferior peribrmance when the channel parameters are known exactiy. 
However, under practical conditions of channel mismatch, the second embodiment often performs better. In terms of 
adaptive implementation, tiiis embodiment requires less information (it does not require explicit channel estimates) but 

IS may be slower to adapt As with ttie first embodiment, the front end of the second embodiment consists again of a bank 
of KL matched f Qters for each of ttie P antennas followed by weighting and combining. However, unlike ttie first embod- 
iment, the order of the real operator and linear combiner(K-by-/C matrix multiplication) are exchanged. In doing so. ttie 
weighting, combining, and linear combiner can be represented by a single complex K-by-KLP mat'ix Wq, The real parts 
of the components in the resulting X-vector are passed to either a decision device or a decoder for further processing. 

20 In ttie adaptive implementation, adaptive aigorittims can be used to ot>tain ttie matrix W^; hence channel estimates are 
not explicitiy required. 

[0008] For quadrature PAM (QAM) data signals, real operators are not required in the detector since the signal con- 
stellation is two-dimensional. In ttils case, a variation on ttie second embodiment would be used which does not use ttie 
real operators but follows ttie liner combiner directly witti tiie appropriate decision device or decoder. 

25 [0009] While the MMSE detection techniques are powerful in and of themselves, their performance can be further 
enhanced by using them witti other multiuser detection techniques such as interference cancellation where interfer- 
ence, instead of being projected away as is done witti linear multiuser detectors, is explicitly subtracted from ttie 
received signal. Interference cancellation can occur prior to or after the MMSE linear con^iner. 
[001 0] Both embodiments of ttie MMSE detector can be generalized for practical systems which utilize, in addition to 

30 the data bearing channels, supplemental channels which may serve as pilot signals. Also, ttie MMSE detectors can 
operate In systems where signals are transmitted with different spreading fectors. 

[001 1 ] A CDMA signal detector in accordance with the invention can provide signif icantiy improved performance rel- 
ative to a conventional signal detector which does not use multiuser detection or which does use multiuser detection 
but does not adaptively mitigate unknown interference. 

35 

Brief Desert ption of the Drawings 

[001 2] The foregoing and ottier features of the invention may become more apparent when the ensuing description 
Is read together witti the drawings, in which: 

40 

Rgure 1 shows a prior art single-user space-time rake receiver enrploying a multi-element antenna array, a bank of 
A/-chip filters at each antenna, and a channel weighter and combiner; 
Figure 2 shows the generic eni)odiment of the invention; 

Rgure 3 shows a first embodiment (Detector A) of the multiuser space-time MMSE detector, 
45 Rgure 4 shows a second embodiment (Detector B) of ttie multiuser space-time MMSE detector, 
Rgure 5 shows an adaptive inplementation of tiie second embodiment; 

Figure 6 shows interference cancellation located prior to ttie linear combiner of ttie first enixxliment; 
Rgure 7 shows interference cancellation located prior to ttie linear combiner of ttie second embodiment; and 
Rgure 8 shows interference cancellation following the linear combiner of the multiuser space-time MMSE detector 
so in accordance witti the invention. 

General Description 

[0013] Consider a K user system where ttie Mh user {k^l ... K) modulates its data sequence b^Q with an N-ch^ 
55 spreading sequence S/^Q. For Detector A, ttie data sequence is pulse amplitude modulated (PAM); for Detector B, ttie 
data sequence can be eittier PAM or quadrature PAM (QAM).. The ti-ansmitted signal undergoes frequency-selective 
fading in ttie channel and arrives at ttie receiver as L time-resolval^e multipath components with complex fading chan- 
nel coefficients c^/(0 ..c^t(0- Assume that ttie receiver is a P-element linear array. If each resolvable multpatti com- 
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ponent arrives as a planar wavefront with angle 6 i^t) with respect to a linear array and if the array spacing is 
suffidently dose (e.g., }J2) so that there is perfect correlation among array elements for a given wavefront. then the 
phase offset of the pth element, with respect to the first is /'j^p(0=exp(n(p-1)sine ^;(0) • The received signal at the 
pth antenna for given synrbd period (and ignoring intersymbo! interference) is: 

K L 
M /.I 

where Ai^ is the amplitude for user k, xf^i is the delay for the Mh user's hh multipath, and np{t) is the additive Gaussian 
noise process which accounts for out-of-cell interference and t>ackground noise. We now make the following assump- 
tions to simplify the analysis. 

(a) The signals are received bit synchronously. 

(b) The time spread is small compared to the symbol period so that intersymbol interference can be ignored.; and 

(c) The phase offsets and channel coeffidents are constant over a symbol period. 

The first assunption will be dropped later. Under these assumptions, the dhip matched-filter output for the received sig- 
nal in equation (1) is a complex A/-vector 

K L 



where S/^^/ is the chip matched-filter N-vector corresponding to S^^f -x;^/) and Hp, is the complex N-vector corresponding 
to the Gaussian noise. The spreading codes are normalized to have unit energy: 



The spreading codes are assumed to be random. However, adaptive implementations of these MMSE detectors require 
the use of short spreading codes which repeat after a few symbol periods. The following notation will be used hereinaf- 
ter: 

^Ki-\\i,V'^Ki.P^^ conrplex P-vector of array coeffidents 

H=[hi •hi,/....hK; v..hK;J complex PxKL array matrix 

vP = c//gg(hi i...hi j^...hK,i...h^ J complex KLP x KL an-ay matrix 

=[Qc,i ■• Or, J complex L-vector of channel coefficis 

C= c//ap(ci...C/c) conplex KLxK channel matrix 

R complex KLxKLconrelation matrix defined by: 



A a diag{A^..Af^) real KxK amplitude matrix 

b^lb^ ,„btfy real AC-vector of data bits 

lu uxu identity matrix 

1 u u-vector of ones 
O^xu X a matrix of zeros 

Ou u-vector of zeros 



The corresponding estimated values for the an^ay and channel coeffidents will be represented with the symbol over 
the symbol for the value. The noise vector is a zero-mean complex Gaussian vector whose distribution can be written 
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in terms of its (component-wise) real and imaginary components: 



5 



X»)J~AkJ 




where we define the real and imaginary operations for matrices and vectors as 



10 



Re(J0=(X+X:)/2 and lh(JQ=(X-X:)/2 



75 where * denotes the complex conjugate. Hence Re(n) and lm(n) are zero-mean Gaussian random vectors whose com- 
ponents have variance and are mutually uncorrelated. 

Conventional (Space-Time Rake) Receiver, (Rgure 1) 

20 [0014] In the context of this array-muKipath channel, the prior art single user detector is a conrelator matched to the 
desired user's conposite array-multipath-spreading code signal. This detector does not account for the presence of 
interferers; however, it is the maximum-likelihood detector if there are no interferers or if they are orthogonal to the 
desired user in anray-code space. As shown in Rgure 1, the detector consists of a bank of correlators 10 at each 
antenna, matched to the KL multipath spreading codes Si j.-.s^ .l -^/c.i -S/c,/.* notation ( S;^/;)in the left-most boxes 

25 of Rgure 1 indicates taking the inner product between Sf^i (the Mh user's spreading code corresponding to the Ah mul- 
tipath delay) and the input vector to the box. The dot In the notation (S/(.^')represents the input. The timing estimates 
are obtained using a separate timing estimation algorithm. The output of the upper left-hand box (which corresponds to 
the first component of the vector z^) is 



where the H superscript signifies the Henmitian transpose (take the conrplex conjugate of each element, then take the 
35 transpose of the resulting vector) of a complex vector. The matched filter outputs for at least a subset of the multiple 
signals, a plurality of the multipath components of these signals, and a plurality of receiver antennas are indicated in the 
drawing at (z^ ... Zp) 

[0015] At reference number 11 , the conrelator outputs are weighted by the complex conjugate of their con-esponding 
channel estimate. More specif k»Ily, the correlator output at the pth antenna for the Ah multipath of user k is weighted 

40 by the complex conjugate of the estimate of the conresponding channel (anrayAmiKipath) coefficient h i^j^c f^j. These 
estimates are obtained using a separate channel estimation algorittim. At each antenna, the L compOTients for the Mh 
user are added (1 2), and tiien tiie resulting P components for the Mt\ user are added 1 3;^ (/c » 1 ... K). Each component 
is then passed to either a decision des^ice or a decoder. For PAM data signals, the decision device for the Mh user (14;^, 
k= 1 ... K) gives as its output the symbol which is closest in the EudUean distance sense to the real part of a Mh vector 

45 component. For QAM data signals, the real operator is not required. 
[001 6] From Rgure 1 . the KLP-vector 



30 




SO 




at ttie output of the matched f Qters can be written as: 



55 
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where we deTine the KLP x KLP matrix 

and the KLP x KL matrix . 9 denotes the Kronecker product operation between two matrices, and the complex noise 
vector is a conplex Gaussian random vector with distrOxition 



10 



Lla<-')J \ki i-In^J 



75 



[0017] Using the feet that 



20 



(ft''e)''SH''c=e*[(ft''H).R]c 



25 (where • denotes the component-wise product between two same-sized matrices), the K vector at the output of the 
channel combiners and the input to the decision devices 14^-1 4k can be concisely written as: 



30 



(3) 



where H denotes the complex conjugate transpose, represents the component-wise product of two same-sized matri- 
ces, and 



35 



MsRejc^KHUjoRjc}. 



40 



The noise vector is purely real and has distribution: 



45 



where 



50 



55 



Ml 



C»[(ft''ft).R]< 



If the data is BPSK modulated, the bit decision for user k is simply the hard limit of the kth component of 
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My«.}:4-5gp(My-i})- 



5 



The corresponding bit error rate Pj^'^^ of the conventional rake receiver is: 



10 




(4) 



IS 



where X^;^;^ is the (k^l^th element of matrix X. The performance comparisons among the various detectors are given in 
terms of BPSK modulated data. t)ut they can also be done for any QAM modulated data in general. 

20 Generic Space-Time Linear iUluitluser Ralce Receiver (Figure 2) 

[0018] The prior art space-time rake receiver described above Is a single user receiver in the sense that for a given 
user, the demodulation uses information only from that user. Since it does not account for the presence of interference 
from other users, its performance suffers. Figure 2 shows a generic multiuser space-time rake receiver. It accounts for 

25 the multiaccess interference using a linear combiner. Two versions of this detector are given. The first, known as detec- 
tor A, uses the real operators before the linear combiner. In this case, reference 20 In Rgure 2 consists of a bank of real 
operators followed by a linear combiner, and each of the components 14^.. .14^ is a slicer which determines the closest 
estimated symbol to its input. The second, detector B. uses a linear combiner represented by a complex matrix followed 
by the real operators. In this case, reference 20 consists only of the linear combiner, and the decision devices consists 

30 of a real operator and slicer. Despite the similarities, the fact that the linear combiners for detectors A and B results in 
different performances and adaptive inplementations. We now describe the two detector embodiments nrore thor- 
oughly. 

Space-Time Linear iVIuitiuser Detector A (Rgure 3) 

35 

[001 9] A set of sufficient statistics for the K users is given in equation (3) as the K-vector Re{yconi^}. The objective of 
the linear MMSE detector is to apply a /Cx K linear combiner to this vector such that the mean-squared en^or between 
its resulting vector and the data vector b is minimized. In other words find the real KxK matrix W4 such that 



This detector in essence strikes the proper balance in coni^ating both residual Gaussian noise, manifested through 
45 n^"^ . and multiple access interference, manifested through the off-diagonal terms of 



40 






50 



It can be shown that the solution is: 



55 
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(6) 



= AM(MA^M''+ct^M)"' 



[0020] Rgure 3 shows a block diagram of detector-A in accordance with the invention. R contains P antennas and a 
bank of matched filters 1 0. The box 3 1 represents a linear combiner defined by the K-by-/C matrix given by equation 
(6). The matrix can be put into words as follows: rt is the product of 3 terms: (a) the diagonal amplitude matrix (A); (b) 
the equivalent array/channel/code correlation matrix (M); and the matrix inverse of the matrix 



which consists of tiie product of the arr^/channel/code conrelation matrix (IV!), the square of the anplitude matrix (A^), 
and the Hemiitian conjugate of the anray/channel/oode conrelation matrix (M^ sunruned with the estimated anray/chan- 
nel/code correlation matrix 



which is weighted by the background noise variance (a^). This matrix minimizes the mean-squared en^or criterion given 
in equation (5), and because of the convexity of the mean-squared enror expression in (5), the MMSE solution (6) rep- 
resents a global minimum. A bank of decision slicers (32^ - 32k) follows tiie linear combiner. 
[0021 ] An attractive feature of the MMSE detector is tiiat it can be inplemented adaptively using well-known adaptive 
algorithms (33) such as least-mean-squares or recursive least-squares. Using a training signal (for exanple, tiie user 
data bits), tiiese techniques can be used to adaptively obtain Wy(. Adaptive implementation is an option if direct calcu- 
lation of W>) is deemed too complex. Note tiiat tiie adaptive implementation requires knowledge of the AC users' spread- 
ing codes, multipath delays, and channel parameters. 
[0022] Assuming that the array and channel estimates are exact 



and anotiier attractive feature is that MMSE detector approaches the decorrelating (zero-forcing) detector (up to a sca- 
lar factor) as tiie residual Gaussian noise approaches zero: 



(MA^M'^+CT^M) 




C=QH=H. 



we can rewrite tiie expression for Wy^ in (6) as: 






In ottier words, from (3) we have that 
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^0 and the multiple access interference is forced to be zera In the inrrplementation of the deconrelating detector A (DD-A), 
we replace given by equation (6) with 



15 



We will assume that C^CtH'^HjJflC is positive definite and hence invertiUe; and it follows that Re{C^[(H^H)i)]C} is 
also invertible. 

[0023] Let us now calculate the bit en^or rate of this first MMSE detector (MMSE A). The K-vector at the slicer inputs is: 



20 



25 where 



30 



Hence the bit error rate for the kth user using the MMSE-A detector is: 



35 



40 



Letting 



45 



(w,MAb)t 



= a-'(rc{c*[(h*h)or]c})"' 



(8) 



and assuming perfect channel estimateSp the bit enror rate for DD-A, is: 



50 



55 
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/r-v)=Gi 



(9) 



10 



[0024] One drawback of this MMSE detector is that the estimation of the channel and an^ay coefficients cannot be 
incorporated into the adaptive algorithm for obtaining W>). The estimates must be explicitly obtained using some other 
IS means (like a training or pilot signal). This leads us to asK Is it possible to design an adaptive MMSE which does not 
require prior knowledge of the array and channel parameters?" The answer Is a"yes." as we will see in the following sec- 
tion. 



20 



45 



55 



Space-Time Linear l\/iultiuser Detector B (Rgure 4) 



[0025] The goal of tiiis subsection is to derive an MMSE detector which can be inplemented adaptively and which 
does not require explicit an^ay and channel estimates. (It will, however, still require timing estimates for the multipatti 
delays.) Figure 4 shows the basic structure of Detector B which exchanges the places of the real operators and linear 
combiner of Detector A. Following tiiis exchange, the resulting cascade of the array combining, multipath combining, 
25 and linear combiner can be collected into a single complex Kx KLP matrix multiplication called Wg. as shown in Figure 
5. Specifically, we define this matrix using the MMSE criterion as: 

W,=argimBj£(|Wz-bfj. 

and it follows from (4) that 

w,=£K)[£(«-)]" 

=AG''R(RGA*G''R+<r'R]'\ (10) 



where 

GsifC 



is tfie ray/channel matrix. Rgure 5 shows a block diagram of ttie MMSE-B detector in accordance with the invention. It 
so contains P antennas and a bank of matched filters 10. The box 50 represents a linear combiner defined by the K-by- 
KLP matrix given by equation (10). This matrix minimizes the mean-squared error critoion given above by 



W, = argimDj£[|Wi-bfJ. 



The last line of equation (10) states that the matrix is the product of 3 terms: (a) the diagonal amplitude matrix (A); (b) 
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the array/channel combiner matrix H^C which is equivalent to the array/channel conft)iner in the conventional detector 
but using the actual parameters and not their estimates; and (c) the matrix inverse of the matrix 

The matrix inverse consists of the product of the correlation matrix 

(R) 



defined following equation (2). the ray/channel combiner matrix G, the square of the amplitude matrix (A^). arKi the Her- 
mltian conjugate of the array/channel combiner matrix summed with the KLP-by-KLP identity matrix weighted by the 
background noise variance (o^. K the data is PAM., a bank of real operators (30<| - 30^) and decision slicers (32^ • 32k) 
follows the linear conri>iner. Othenwise, if the data is QAM, only the complex decision slicers (14^ -14k) required. For 
PAM data, Detectors A and B, while similar In structure, are not equivalent due to the different placement of the real 
operator. 

[0026] To calculate the bit enror rate of Detector B for BPSK data, the /C-vector at the slicer inputs is 

Ke[ya}=R<W^=F£(W^C)Ab+n|, 



where 

Hence the bit en'or rate for the Mt\ user using the MMSE-B detector is: 




[0027] By concatenating the P received AZ-vectors into a single NP-vecXor 

we can combine the operations of the matched filter bank with the linear conAiner W0 to create a single K-by-A/P 
matrix which multipltes r. As shown in Rgure 6, the detector structure for the Mh user can be reduced to the Mh row of 
this K-by-NP matrix 60 followed by the real operator 62 and decision slicer 63. Hence the MMSE-B detector can be 
implemented with a simple tap-weight filter archrtecture. Hs conresponding deconrelating detector (described below), the 
MMSE-A detector and the DD-A detector can likewise be implemented with this architecture. Like the MMSE-A detec- 
tor, this version of the MMSE-B detector can be inplemented adaptively using well-known adaptive algorithms (61). 
However, MMSE-B has the advantage of not requiring the array and channel estimates. The tradeoff is that the adap- 
tation may be slower because there are more taps to adjust 

[0028] While direct analytic conparisons between the two detectors is difficult, we will show that under perfect chan- 
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nel and array estimates, the performance of MMSE-A is uniformly superior to that of MMSE-B as the noise floor 
approaches zero. Rrst we show that the corresponding decorrelating detector for MMSE-B, DD-B, simply replaces 
with 

A-f(H^cfRH''Cr(H^cf. 

(Recall that we have already made the assumption that 

10 

(H''cfte^C=C^[(H''H)oR]c 

IS 

Is invertible.) 
Proposition 1: 
20 [0029] 

lim A(H^C)''[RH^CA'(H''cf + (t'I^]"' j - A-|(H^cf RH^c|"(H^cf 



25 



Proof: Multiplying both sides by 

30 



RH''CA'(rc)%CT^I^, 



35 we have that 



40 



= A(H''cf + lm(<7' A-|(H''cf fe^c|''(H''C)'') 
«A(H''cf. 



45 

Hence, from the fact that 

so 



(rcffiH^c=c"[(rH)oR]c 



and the fact that the array and multipath combiners in MMSE-A can be represented by (H^D)^, the DD-B is the bank 
of matched fitters, followed by the combiners, followed by the deconrelator (C^[(H^H)A]C)'\ followed by the real oper- 
55 ators. and followed by the slicers. The DD-A is nearly Identical except that the order of the decorrelator and the real 
operators are switched, and the deconelator is (Re{C^[(H^H).R]C})'^ 
By r^laclng W3 with 
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5 



in (1 1), the bit error rate of DD-B is: 



10 



4 




(12) 



IS 



20 [0030] We would now like to show that DD-A has a uniformly lower bit error rate than DD-B. 
Proposition 2: P^^^ (o) ^ P j?^^ (a) 

[0031] Proof (for which the inventors wish to thank Emre Telatar of Lucent Technologies* Math Sciences Center for 
25 his timely and elegant insight): 
[0032] Define 



which is Hermitian and assumed to be positive definite. It follows that [Re(RI)]'^ is also positive definite, and it is suffi- 
cient to show that the kth diagonal element of [Re(IBI)]'^ is greater than or equal to the kth diagonal element oiW,Jllt' 
V (Note that since KR is Hermitian. its inverse is also Hermrtian, hence the diagonal elements of the inverse are real.) 
35 Chose y to be real and set 



30 





40 



and note that X is real. Thus, 



x^Mjl = x^Re(M)x = y[Re(M)] 



4S 



and also 



so 




Using a generalization of Bergstrom's inequality which states that for any real x, real y, and positive definite ffl, 



55 




we have that 
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y'[R«(M)]"'yy^(M)'y i (/[ReHf'')' 
y'(5!)-'yiy'(Re(5i)fy. 

10 

and the result follows by using the kth unit vector for y. Despite this relationship, it turns out that the performance of DD- 
A and DD-B are actually quite similar. More Importantly, it has been shown that the performance of MMSE*B can be 
better than MMSE-A for the practical case of inexact anray and channel estimates. 

15 Extensions of MMSE Detectors to the Asynchronous MuWrate Case 

[0033] So far, we have assumed bit synchronous operation of these linear multiuser detectors. We now show how 
their operation can be extended to the bit asynchronous case. First, we consider the operation of a decorrelating detec- 
ts in a synchronous single-rate environment. Similar reasoning applies to MMSE detectors; however, for pedagogical 
20 purposes, we focus on the decorrelating detector. Suppose that there are two users with respective spreading codes 
Si(0 and S2(t). Table A below shows the relationship between the two codes in time, with the thin vertical lines repre- 
senting symbol boundaries. The received signal is r(t)^A^s^{f)b^ + ^252(0^2 + ^^(0 where Aff Is the kth user's 
amplitude, bf^ Is the Mh user's data bit, and n(0 is the additive Gaussian noise. 



• • 




• • 


• • 




• • 



Table -A: 2 user synchronous single-rate system 
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',(0 


• • • 


• 




s,it-T) . . . 


»-B: 2 user asynchronous single-rate s^tem 


• • 




• • • 


• • 


',(0 S,(t'T) 


Ji(/-27) . . . 
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Table • C: 2 user asyndironous multi-rate system 



25 



30 



[0034] If user 1 is the desired user. Its decorrelating detector for a given symbol interval is the matched filter for code 
Si(0 projected into the null-space of S2(Q. For the asynchronous single-rate case, shown in Tat)te - B, the decorrelating 
detector for user 1 is projected into the null-space spanned by the linear combinations of user 2's codes which overlap 
during that symbol Interval. If we let s ^(0 be the portion of $2(0 which overlaps with s^{t) and let (f- 7) be the portion 
of Szit'J) which overlaps with s<|(0, then, assuming BPSKdata modulation, the null space Is spanned by 

Sziths^it'T) andsf(0-S2(M). 

This idea can be extended to the asynchronous mutti-rate case shown in Table - C. Using similar function definitions, 
the decorrelator for user 1 lies in the null space spanned by 



35 



s^(0-S2(f-T)+4(^-2r). and s^{tys^{hT)'S^{h2n 



As the disparity between data rates increases, the effectiveness of a linear multiuser detector decreases since it will be 
40 constrained to a more restrictive subspace. The reasoning above can be applied to MMSE detectors by similarly com- 
bining the appropriate subspaces of the Interfering waveforms. 

MMSE Detectors and Interference Cancellation 

45 [0035] TTie MMSE techniques described here can be used inconjunction with a non-linear multiuser detection tech- 
nkiue known as interference cancellation. Whereas the linear multiuser detectors rely on subspace projections to miti- 
gate interference, the interference cancellers directly subtract off the interference. Two types of interference cancellation 
can be used: pre-combiner. where the cancellation occurs before the linear combine, and post-combiner, where it 
occurs after. Pre-combiner cancellation is used for canceling interference which is not accounted for by the linear com- 

50 biner. For example, we may choose to cancel a high-powered, high-data rate signal from within the cell instead of 
accounting for it with the linear combiner. Post-combiner cancellation on the other hand Is used for refining the symbol 
estimates made by the linear combiner. Using preliminary syn*ol estimates from the linear multiuser detector, the sig- 
nals are reconstructed and subtracted from the received signals to form enhanced signals which are used for a second 
stage of syntK>l estimation. In both interference cancellation techniques, capacity can be potentially Increased by sup- 

55 pressing the interference. 

[0036] Rgure 7 shows a block diagram of a pre-combiner interference canceller in accordance witii Detector A shown 
in Rgure 3. At the Mh Input to the linear combiner W^, there Is interference due to users y « 1 ...K, j ^ k, which will be 
accounted for with Wy^. However, there may also be othenwise unaccounted Interference, for exanple, high-powered 
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interference from within the cell. This Interference could be subtracted from the input to Wy^ if the spreading codes, mul- 
tipath delays, data bits, arKi channel parameters at all antennas are known. In this case, the baseband signal for the 
interference could be reconstructed and its contribution to the Mh input of W>) could be estimated (by passing the recon- 
structed signal through the processing chain for user k) and subtracted at the point designated 70/^ (/c = 1 ... K) Contri- 

5 buttons from other interference for this or other inputs can be similarly computed and subtracted. 

[0037] Rgure 8 shows a block diagram of a pre-combiner interference canceller in accordance with Detector B from 
Rgure 4. Using the same method as described above for the Detector A pre-combiner interference canceller, the inter- 
ference which is not accounted for by can be subtracted from each of its inputs at the point designated 80. 
[0038] Rgure 9 shows a block cfiagram of a post-combiner interference cancellation receiver in accordance with the 

10 invention. This receiver irK^ludes P antennas and has three stages, designated 90. 91 arxi 92. The first stage 90 
receives baseband signals fp, and makes preliminary symbol estimates for all /C users with either a MMSE-A detector 
of Rgure 1 or a MMSE-B detector of Rgure 2. The second stage 91 uses the preliminary symbol estimates from the 
first stage 90 and knowledge of the K users' spreading codes, delays, and channel parameters to reconstruct the base- 
band received signals for each of the users. The second stage 91 subtracts the multiaccess interference with respect 

IS to a desired user k from the received signals .. .fp to f onm the following enhanced received signals for user k. 



20 



K L 



[0039] Note that the multipath interference with respect to user k can also be removed. We will not describe this option 
in detail because its relative gains will generally be negligible unless the number of users K is very small. Ideally, if the 
25 preliminary symbol estimates and the channel estin^ates are perfect, there will be no multiaccess Interference in the 
enhanced received signals. However, this will generally not be the case in practice. In the third stage 82, tiie enhance 
received signals are processed witii a conventional space- time receiver shown in Rgure 1 to generate flie final symbol 
estimates b J/^) for the K users. 

[0040] The post-combiner interference cancellation procedure described above can be continued for mutttpie itera- 
30 tions by repeating stages 91 and 92 indefinitely Performance gains due to repeated iterations diminish after the first two 
or three iterations. To reduce the complexity, interference cancellation of botti types can be performed on a subset of 
users. 

[0041] What has been described is deemed to be iflustrative of the principles of ttie invention. In particular, it should 
be noted that the apparatus and methods of the invention may be implemented tiirough various technologies, for exam- 
35 pie, tiirough ttie use of large scale integrated circuitry, application specific integrated circuits and/or stored program 
general purpose or special purpose computers or microprocessors using any of a variety of computer-readable media. 

Claims 

40 1 . A method for detecting the data signals modulated respectively by multiple direct-sequence spread spectrum sig- 
nals, comprising the steps of 

a) generating matched filter outputs for at least a subset of the multiple signals, a plurality of the multipath com- 
ponents of these signals, and a plurality of receiver antennas; 
45 b) developing coherent channel estimates; 

c) weighting and combining the matched filter outputs using ttie coherent channel estimates; and 

d) estimating the respective data symbols for at least a subset of ttie spread spectrum signals while siq^press- 
ing multiple access Interference. 

so 2. The method of claim 1 wherein said coherent channel estimates are developed from said matched filter outputs. 

3. The mettiod of daim 1 wherein said coherent channel estimates are developed from an auxiliary signal. 

4. The method of daim 1 wherein multiple access interference Is suppressed with a linear combiner and wherein a 
55 decision device optimally estimates ttie respective data syntels given the linear combiner outputs. 

5. The method of claim 4 wherdn ttie steps of generating matched filter output; of combining the matched filter out- 
puts, and of utilizing a linear combiner are implemented using a modified tap-weight filter. 
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6. The method of claim 4 wherein the linear combiner employs a minimum mean-squared error criterion. 

7. The method of daim 6 wherein the minimum mean-squared error criterion designates a linear combiner for the Mh 
user which minimizes the expected squared enror between the conf4)iner output and the Mh user data symbol. 

8. The method of claim 7 wherein the linear combiner forms a matrix Wa from the complex block-Toeplitz conrelation 
matrix 

a. 

the diagonal amplitude matrix A, the conplex ray/channel matrix G, and the background noise variance o^. 

9. The method of claim 7 wherein the linear con*iner forms a matrix Wg from the complex block-Toeplitz congelation 
matrix 

the diagonal amplitude matrix A, an estimate of the complex anray/channel matrix 

G. 

and the background noise variance o^. 

10. The method of claim 4 wherein the linear combiner emf^oys a zero-forcing error criterion. 

11. The method of claim 10 wherein the linear combiner fbnns a matrix from the complex block-Toeplitz congelation 
matrix 

the complex ray/channel matrix G, and the diagonal amplitude matrix A. 

1 2. The method of daim 1 0 wherein the linear combiner fbmns a matrix from the complex block-Toeplitz conrelation 
matrix 

8. 

an estimate of the complex anray/channel matrix 

G, 

and the diagonal ampfitude matrix A. 

13. The method of daim 4 wherein the data signals are pulse amplitude modulated, and the real components of each 
of the combined matched filter outputs are extracted prior to the linear combiner. 

14. The method of claim 13 wherein the linear combiner employs a minimum mean-squared enror criterion. 

15. The method of claim 14 wherein the minimum mean-squared error criterion designates a linear combiner for the 
kHh user which mininvzes the expected squared enror between the combiner output and the Mh user data symbol. 
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16. The method of daim 1 5 wha'ein the linear combiner fomns a matrix Wy^ from the real equivalent anray/channel/oode 
correlation matrix Re(M). the diagonal amplitude matrix A, and the background noise variance o^. 

17. The method of daim 15 wherein the linear coni)iner fomis a matrix W>) from the real part of an estimate of the 
equivalent array/channel/code correlation matrix 



the diagonal amplitude matrix A, and the background noise variance o^. 

18. The method of daim 13 wherein the linear conrtbiner employs a zero-fordng error criterion. 

1 9. The method of claim 1 8 wherein the linear combiner zero-fordng criterion indudes fonming a matrix from matrix 
inverse of the real equivalent array/channel/code congelation matrix Re(R/0. 

20. The method of daim 18 wherein the linear combiner zero-fordng criterion indudes forming a matrix W^^ from the 
matrix inverse of the real part of an estimate of the equivalent an^ay/channel/code coaelation matrix 



21. The method of daim 1 wherein the direct sequence spread-spectrum signals are spread with different spreading 
gains. 

22. The method of claim 4 wherein, prior to employing the linear combiner, interference cancellation is performed on at 
least a subset of the combined matched filter outputs. 

23. The method of daim 4 wherein the linear combiner is followed by interference cancellation on at least a subset of 
the linear combiner outputs. 

24. Apparatus for use in detecting at a receiver in a wireless system with multiple transmitted direct-sequence spread 
spectrum signals, comprising: 

a) means for generating matched filter outputs for at least a subset of the multiple signals, a plurality of the mul- 
tipath components of these signals, and a plurality of receiver antennas; 

b) means for developing coherent channel estimates from the matched filter outputs; 

c) means for weighting and confining the matched filter outputs using the coherent channel estimates; and 
cO means for estimating the respective data signals for at least a subset of the spread spectrum signals while 
suppressing multiple access interference. 

25. A receiver comprising: 

a) space-time rake receiver having: 



i) a bank of correlators for resolving K code-spread data signals Si,S2,^„,Sf^ bearing respective data 
b'i,b2,...,bfc each with at most L resolvable multipath components with conplex amplitude Cj^ {k = 

1.2 K; / s 1,2,...,L), arriving at a P element antenna array with complex amplitude /7^/^ (p « 1,2,...,P) 

into at most KLP correlator outputs corresponding to at most L multipath signals for each of the AC users 
at each of the P antenna array elements; 

ii) means for weighting the output of each correlator accoiding to the conplex conjugate of an estimate of 
the coefficient for a respedive channel and path {h )^ ^ ^ c \ f): 

iii) means for summing together the weighted resulting LP complex values for each user, and 

iv) means for forming a K-vedor of the AC resulting sununalions; with 
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b) means for applying a linear coiTt}iner to the /C-vector to suppress multiple access interference; and 

c) means for processing the resulting signals to produce estimates of the respective data symbols. 

26. A receiver according to daim 25 wherein said means for the linear combiner comprises means for forming the prod- 
uct of: 

the diagonal amplitude matrix (A); 

the equivalent anray/channel/code correlation matrix (M); and 

the matrix inverse of the matrix 



where is the hermitian conjugate of the anray/channel/code correlation matrix, 

(M) 



is the estimated array/channel/code correlation matrix and (o^) is the background noise variance. 

27. A receiver according to claim 25 wherein the real parts of the K resulting summations are extracted prior to the lin- 
ear combiner. 

28. A receiver comprising: 

a bank of correlators for despreading K code-spread data signals Si.S2....*Sk bearing respective data 
bi,/>2. 3dch with at most L resolvable multipath components to form a /C-vector, 
means for applying a linear combiner to the K-vector to suppress multiple access interference; and 
means for processing the output of the linear combiner to produce estimates of the respective data symbols. 

29. A receiver according to claim 28 wherein said means for the linear combiner comprises forming the product of: 

the diagonal amplitude matrix (A); 

the hermitian transpose of the anray/channel matrix; and the matrix inverse of tiie matrix 



where 

s 



is ttie complex block-Toeplitz conrelation matrix, is the background noise variance, and I^lp is the identity 
matrix whose dimension is KLP x KLP. 

30. A method of demodulating K direct sequence code-spread data signals sutqected to L-pati) frequency-selective 
bding, received at an anray of P antennas comprising: 

despreading with a bemk of filters notched to L multipath delayed replicas of each of the K signals received at 
each of the P antennas; 

weighting the despread signals with the complex conjugate of an estimate of the ray/channel coeffident; 
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summing together the weighted, despread signals to form a K-vector; 

applying a linear transfonn to the K-vector to sif)pres8 multiple access interference; and 

processing the resulting signals to produce estimates of the respective data 8ynri)ols. 

31 . A method according to daim 30. wherein said linear transform is the product of: 

the diagonal amplitude matrix (A); 

the equivalent anray/channel/code congelation matrix (M); and 

the matrix inverse of the matrix 

(MAW+(T^M), 



where h/h is the hermitian conjugate of the anray/channel/code correlation matrix, 

H 

is the estimated anray/channel/code conrelation matrix and (o^) is the background noise variance. 

32. A method of demodulating K direct sequence code-spread data signals subjected to L-path frequency-selective 
fading, received at an ray of P antennas comprising: 

a) despreading with a bank of filters matched to L multipath delayed replicas of each of the /C signals received 
at each of the P antennas to produce a /C-vector, 

b) applying a linear transform to the K-vector to combine the ntatched filter outputs and suppress multiple 
access interference; and 

c) processing the resultant signals to produce estimates of the respective data symbols. 

33. A method according to daim 32 wherein said means for the linear conlblner comprises forming the product of: 

d) the diagonal amplitude matrix (A); 

e) the hermitian transpose of the anray/channel matrix (G'^; and 

f) the matrix inverse of the matrix 

where 

R 

is the complex block-Toepiitz con^elation matrix, o^ is the background noise variance, and Iklp is the identity 
matrix whose dimension is KLP x KLP. 

34. A recaver for use with direct-sequence spread spectrum signals, con^rislng: 

a rake receiver with array processing; and 

a multiuser detector receiving an output of said rake receiver. 

35. A receiver according to daim 34 wherein said rake reiver with array processing indudes: 

a bank of f ilters matched to ttie multipath timing delays of the multiple spread spectrum signals; and 
a drcuit for weighting the matched filter outputs with the conresponding anray/channel estimates. 
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36. A system for detecting the data signals nxxlulated respectively by nultipie direct-sequence spread spectrum sig- 
nals, comprising: 

a plurality of space diversity antennas; 

a plurality of matched filters associated with said antennas for provkiing a plurality of outputs (z^ ... Zp) for at 
least a sut)set of the multiple signals, a pluraTrty of the muttipath components of these signals, and a plurality 
of receiver antennas; 

a circuit for obtaining coherent channel estimates {hf^jj^c f^J) from the matched fOter outputs; 

a second drcut for weighting and combining the matched filter outputs using the coherent channel estimates; 

and 

a third circuit for estimating the respective data symbols for at least a subset of the spread spectrum signals 
while suppressing multiple access Interference. 

37. A receiver comprising a processor for executing software, said software being in a computer readable medium and 
being for detecting the data signals modulated respectively by multiple direct-sequence spread spectrum signals, 
said software causing said processor to perform at least the steps of: 

a) generating matched filter outputs for at least a subset of the multiple signals, a plurality of the multipath com- 
ponents of these signals, and a plurality of receiver antennas; 

b) developing coherent channel estimates from the matched filter outputs; 

c) weighting and confining the matched filter outputs using the coherent channel estimates; and 

(0 estimating the respective data symbols for at least a subset of ttie spread spectrum signals while suppress- 
ing multiple access interference. 



EP 0 978 951 A2 




EP0978 951 A2 




EP 0 978 951 A2 




EP 0 978 951 A2 




EP 0978 951 A2 



FIG. 5 
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